Extending the Velocity-dependent One-scale String Evolution Model 
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We provide a general overview of the velocity-dependent one-scale model for cosmic string evo- 
lution and discuss two further extensions to it. We introduce and justify a new ansatz for the 
momentum parameter k, and also incorporate the effect of radiation backreaction. We thus discuss 
the evolution of the basic large-scale features of cosmic string networks in all relevant cosmological 
scenarios, concentrating in particular on the 'scaling' solutions relevant for each case. In a compan- 
ion paper, we show, by comparing with numerical simulations, that this model provides an accurate 
description of the large-scale features of cosmic string networks. 



I. INTRODUCTION 



o. 

The velocity-dependent one-scale (VOS) model provides the most convenient and reliable method by which to 
calculate the large-scale quantitative properties of a string network in cosmological and other contexts jl] |). It is 
widely used for making quantitative predict ions of the potential observational implications of cosmic strings [0. 
Given its simplicity, it is remarkable how well the VOS model performs when tested against high resolution numerical 
simulations of string networks H . It is well-known that string evolution is a complex physical process with a build-up 
■ of small-scale structure on the strings, which is very computationally demanding to model accurately ■ Analytic 
approaches like the VOS model, abandon the possibility of describing the statistical physics of the string network 
accurately and concentrate instead on its thermodynamics. In other words, a small number of macroscopic quantities 
are selected and the microscopic string equations of motion are used to derive evolution equations for these averaged 
quantities. The price to be paid in this approach is that the averaging process introduces phenomenological parameters 
r"| ■ whose values are not specified by the model itself. Instead, one must still fix these parameters by direct comparison 
CIh! with numerical simulations. 

The VOS model is a generalization of the 'one-scale' model pioneered by Kibble |J (see also ref. |T(|) which describes 
' string motion in terms of a single correlation length L. By incorporating a variable rms velocity v, the VOS model 
extends its validity into early regimes with frictional damping and across the important matter-radiation transition, 
thus giving a quantitative picture of the complete history of a cosmic string network. Other analytic approaches to 
. , string evolution have attempted to incorporate the additional small-scale structure seen in numerical simulations. This 
includes a 'kink-counting' model ipTfl , a functional approach a 'three-scale' model and a 'wiggly' model fli) . 
While these are important for characterising detailed network features, they introduce a significant number of further 
phenomenological parameters which must be fixed by simulations (and which remain rather uncertain). Nevertheless, 
for describing the large-scale properties of a long-string network, the VOS model has proved to be sufficient for a good 
quantitative fit using only a single parameter, the loop chopping efficiency c. 

The purpose of the present paper is, first, to provide a concise exposition of the VOS string evolution model. We 
summarise how it can be applied to describe cosmological string evolution, including late times with a cosmological 
constant, and we present the very different histories of both GUT- and electroweak-scale strings. Secondly, we propose 
an improvement of the VOS model by presenting a new ansatz for the momentum parameter k, which we justify both 
analytically and numerically. Thirdly, we present a further extension incorporating radiation backreaction, which 
provides small corrections to the cosmological scaling laws and which also compares favourably with published results 
of global string simulations. Finally, we review generalizations in a curved FRW spacetime, giving some further 
asymptotic scaling solutions. We report on detailed comparisons between numerical string simulations and the VOS 



model elsewhere [p 15 1 
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II. THE VOS STRING NETWORK MODEL 



The velocity-dependent one-scale model has been described in considerable detail elsewhere llH^Jl7r|l9|l i so here we 
limit ourselves to a brief summary which highlights the features that will be important for what follows. Also for 
simplicity, we will only discuss the evolution of the long string network, even though this formalism is also applicable 
to the loop population, mutatis mutandis. We will discuss this case in detail elsewhere Ea|. 



by 



A. The averaged evolution equations 

Averaged quantities which we could use to describe the string network are its energy E and RMS velocity v defined 



o _ fx 2 eda 



E = paM / eda , v A = J r , , (2.1) 
J J eda 

where the string trajectory x(tr, t) is parametrised by the worldsheet coordinates a and t and the 'energy density' 
e(cr, t) gives the string length per unit a along the string. 

Any string network divides fairly neatly into two distinct populations, long (or 'infinite') strings and small closed 
loops with corresponding quantities denoted by a subscript oo and £ respectively. The long string network is a 
Brownian random walk on large scales and can be characterised by a correlation length L. This can be used to replace 
the energy E^ = PooV in long strings in our averaged description, that is, 

p_ 

Poo - L2 • 

A phenomenological term must then be included to account for the loss of energy from long strings by the production 
of loops, which are much smaller than L. A 'loop chopping efficiency' parameter c is introduced to characterise this 
loop production as 

dpoc \ ~ Poo ro . 

In this approximation, we would expect the loop parameter c to remain constant irrespective of the cosmic regime, 
because it is multiplied by factors which determine the string network self-interaction rate. 

From the microscopic string equations of motion, one can then average to derive the evolution equation for the 
correlation length L, 

2^=2HL(l + v 2 oo ) + jvl + cv 0O , (2.3) 



where H is the Hubble parameter and if is a friction damping length scale. The first term in (2J3) is due to the 
stretching of the network by the Hubble expansion which is modulated by the redshifting of the string velocity. The 
second term is due to frictional interactions by a high density of background particles scattering off the strings. The 
friction length scale if (defined in ref. (!]]) typically depends on the background temperature as If ~ l(l)pT~ 3 , so 
that it grows with the scale factor as a 3 . It usually becomes irrelevant after a time t* ps (Gp)~ 1 t c (with t c being the 
epoch at which the network was formed) which i s a very short time for GUT-scale strings but can be as late as to 



for electroweak strings. Note also that equation (2.3) is valid for an arbitrary flat FRW model with H given by the 
Friedmann equation, 

H 2 = (^j = Hi (Q RQ a- 4 + fi M0 a' 3 ) + ±A . (2.4) 

where fipio and ^mo are the fractional radiation and matter densities today at to, A is the cosmological constant and 
we take a(to) = 1. 

One can also derive an evolution equation for the long string velocity with only a little more than Newton's second 
law 



dVcj , 2 \ 



~k 






{ 2H + Tf) v -_ 


L ~ 





(2.5) 
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where k is called the 'momentum parameter'. The first term is the acceleration due to the curvature of the strings 
and the second damping term is from both the expansion and background friction. The parameter k is defined by 



k _ ((l-x 2 )(x-u)) 
v(l — v 2 ) 



(2.6) 



with x the microscopic string velocity and u a unit vector parallel to the curvature radius vector. In previous work 
p|,p|p^j3]| , we left k as a second phenomenological parameter, while pointing out that it is related to small-scale 
structure and also demonstrating specified asymptotic dependencies on the velocity. In the next section, however, 
we justify an accurate ansatz for k which removes this additional freedom. For most relativistic regimes relevant to 
cosmic strings it is sufficient to define it as follows: 



2V2 1 - 8l 



l + 8v 6 



(2.7) 

2^/2/tt with a more complicated 



In the extrem e fr iction-dominated case (v — > 0), we have the nonrelativistic limit fc„ 
ansatz than (2.7) interpolating between these limits for intermediate regimes. 

Finally, we end this summary by noting that the VOS model has been extensively compared with the results of 
numerical simulations J^PPlJi*! and shown to provide a good fit to the large-scale properties of a string network. 
In particular, it matches well the evolution between asymptotic regimes as a network passes through the matter- 
radiation transition. Comparisons with numerical simulations confirm the constancy of the only free parameter, the 
loop chopping efficiency c, and fix its value to be || 



c = 0.23 ±0.04. 



(2.8) 



The VOS model for any flat FRW cosmology, then, consists of the evolution equations (|2.3|) an d (2.5) with the 



parameters c specified in ( |2.8P and k given by (2.7) (or a more accurate general expression given below) and the scale 
factor a satisfying the Friedmann equation (2.4). 



B. Scale-invariant solutions 



We now start to use the VOS model to provide a general overview of the evolution of string networks in various 
cosmological scenarios. First we analyse some basic late-time properties of cosmic string networks, neglecting the 
effect of friction due to particle scattering. A crucial question is whether or not they can reach a 'scale invariant' 
attractor solution which is required, among other things, for a Harrison-Zel'dovich spectrum of primordial density 



fluctuations to be generated. This can be discussed by analysing the VOS equations ( [2.3] ) and (2.5). 

Scale- invariant solutions of the form L oc t, L oc H^ 1 or L oc dn, together with Voo = const., only appear to exist 
when the scale factor is a power law of the form 

a(t) oc t 13 , (3 = const. , < f3 < 1 . (2.9) 

This condition implies that 

Loctac H^ 1 ocd H , (2.10) 

with the proportionality factors dependent on (3. It is useful to introduce the following useful parameters to describe 
the relative correlation length and densities, defining them respectively as 

L = -yt, ( = 7 ~ 2 = P oot 2 /n . (2.11) 

By looking for stable fixed points in the VOS equations, we can express the actual scaling solutions in the following 
implicit form: 

k(k + c) k{l-f3) 

where k is the constant value of k(v) given by solving the second (implicit) equation for the velocity. Although it 
may not be obvious by inspection, it is easy to verify numerically that this solution is well-behaved and stable for all 
realistic parameter values. 
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If the scale factor is not a power law, then simple scale-invariant solutions like (2.15) do not exist. Physically 
this happens because the network dynamics are unable to adapt rapidly enough to the changes in the background 
cosmology. A prime example of this is, of course, the transition between the radiation and matter-dominated eras. 
The evolution of a GUT-scale network shown in fig. |] illustrates asymptotic regimes in which string evolution 
is scale- invariant, as well as the matter-radiation transition where it is not. Note that for realistic cosmological 
parameter choices, a string network today is still only slowly approaching its asymptotic matter density. Since, the 
cosmological importance of the changes in the network properties during the matter-radiation transition cannot be 
over-emphasised, it is important to calculate these accurately either with direct numerical simulations or with the 
VOS or similar analytic model. The same is also true for late time curvature or cosmological constant domination. 



C. Friction-dominated scaling solutions 



During friction-dominated epochs one has different 'scaling' solutions. However, these are no longer 'scale-invariant', 
since in this case the network retains a memory of its initial conditions, and in particular of the epoch of formation. 
This can be conveniently expressed by a parameter 



tpi 



1/2 



(2.13) 



which is essentially the value of the ratio of the damping terms due to friction and Hubble damping, measured at the 
epoch of string formation. 

Thus in realistic cosmological contexts one can have two different regimes @-||. The first is a 'stretching' regime, 



L 
17 



1/2 



t 



(2.14) 



which is an early-time, transient period which will occur when the initial string density and velocity are sufficiently 
low — for example, as a result of a slow first-order phase transition. In this case the network starts out with a correlation 
length significantly larger than the damping length and so is 'frozen', and is conformally stretched. However the 
damping length is growing as if oc i 3 / 2 , so it quickly catches up with it, ending this regime. However, this can last 
for many orders of magnitude in time for electroweak-scale networks. Although this is not cosmologically relevant 
except for extremely light strings, the analogous regime in the matter-dominated case would be 



L oc t 



2/3 



(2.15) 



The true attractor solution for a friction-dominated epoch, which follows the stretching regime (if this exists) is the 
Kibble regime, which in the radiation era has the form 



L 



30 



1/2 



5/4 



20(c+k nr ) 



1/2 



1/4 



(2.16) 



where k nr is the value of the momentum parameter in the nonrelativistic limit given above. In this case the correlation 
length stays halfway between the damping length and the horizon length. Again there is a matter era analogue, which 
as the form 



L oc t 3/2 , v oc t 1/2 . 



(2.17) 



but this is rarely relevant cosmologically. In the extreme case of an electroweak scale network, friction domination 
ends after radiation-matter equality, but not far enough away from it for this regime to be reached. 

These points are illustrated in fig. || where we see the evolution of an electroweak-scale string network. It is 
interesting to note the scaling behaviour for two cases with very different initial conditions from extreme first- and 
second-order transitions. The high density strings from a second-order transition quickly approach the Kibble scaling 
regime (2.16) disc ussed above. However, the low density strings from a first-order transition begin in a distinct 
stretching regime ( 2.14 ) and persist in it with their density falling slowly until it matches that for the attractor Kibble 
regime. In this case, the string network retains a 'memory' of its initial density for about ten orders of magnitude 
in cosmic time. However, even during the friction-dominated regime the network is able to erase this memory once 
the Kibble regime is reached. On the other hand, the memory of the epoch of formation is not erased — one could 
in principle recover it by measuring the parameter 6. This can also be seen for GUT strings in fig. [j], although the 
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network does not have time to relax into a definite scaling regime before friction-domination ends. For electroweak 
strings, there are also interesting departures from scaling behaviour at the matter-radiation transition and the network 
remains friction-dominated until about three orders of magnitude in time afterwards. Also note that in the Q m = 1 
case the strings only reach the relativistic regime at about the present time, and in the observationally preferred case 
Q rn = 0.2, J7a = 0.8 they are always non-relativistic. This point is crucial, among other things, for a quantitative 
analysis of the evolution of superconducting strings and vortons . 



D. A cosmological constant 

We can also use the VOS model in a flat background to discuss the domination at late times by a cosmological 
constant A [ p2|^9[ , a model for which there appears to be growing observational evidence. In the extreme asymptotic 
case when the universe is inflating we have a oc exp (Ht) with H = y/A/3. The network will 'freeze out' and will 
simply be conformally stretched, that is, 

Lena, Uoo oc a -1 , (2-18) 
where, as soon as the strings become nonrelativistic fc nr = 2\/2/7r, their product satisfies 

L Voo = —H. (2.19) 

Of course, this solution will only apply at early times actually during inflation. At the present time we will only be 
slowly approaching a new stretching regime, so we have to solve the VOS model explicitly. This is shown for a model 
in which f2^ = 0.8, as a dashed line at late times in figs. as well as in detail for GUT-scale strings in fig. |[ It 
is clear that there is a significant fall in the string density and velocity, an effect which, for example, would affect the 
large-angle anisotropies in the cosmic microwave sky. The evolution of the string network is clearly not scale-invariant 
during any period after t eq . 



III. THE MOMENTUM PARAMETER 



Having introduced the VOS model and some of its key cosm ologi cal implications, we now turn to a m ore detailed 



discussion of the so-called 'momentum parameter' k defined in (2.6) and which is important for solving (2.5). String 
velocities are determined by the current acceleration to which they are subjected from the local string curvature, 
as well as their 'bulk' momentum left-over from previous accelerations. Heuristically, we can imagine separating the 
velocity into these the curvature 'c' and bulk 'p' contributions as x = x c +x p . In the extreme friction-dominated limit, 
the velocity is entirely due to the curvature and reaches a limiting average velocity set by the friction length scale, 
v c = (x^) 1 / 2 = lf/L. However, as the velocity increases towards relativistic values we can expect the momentum 
contribution to become larger. Let us suppose that their relative contribution is proportional to some power law of 
the total velocity v, that is, v p /v c oc v a where a is clearly greater than unity. In this case, one can find after some 
straightforward (although tedious) manipulations that the following approximate relation holds for the momentum 
parameter k 

k~ fJL- , (3.1) 

1 + 2 a v 2a v ' 

In this we have also used the fact that flat spacetime analytic calculations have shown that 

fc(l/ v / 2) = 0, (3.2) 



which holds exactly. Note that the above expression means that (3.1) can also be approximately written as 

, 1 - {Vp/Vcf 



1 + (Vp/V c ) 



2 



(3.3) 



We can determine a by studying the well-known ]23| ] helicoidal string solution in flat space but perturbed by a 
frictional force. This solution is 
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x = ( Asincr(cost + rf), ^4coscr(cost + 77), \f \ — A 2 a 



(3.4) 



where < A < 1 and n is a small perturbation, which vanishes if there is no friction. Here, A — 1 corresponds to a 
circular loop, while A = is a static straight string. The evolution equation for the perturbation r\ has the form 



?7 + 2r) 



sin t cos t 
I- A 2 sin 2 t 



(1 - A 2 ) sin 2 t 



cos 2 £ 



1 - A 2 sin 2 t 



sini 



A 2 sin 2 i) 



(3.5) 



where if is the friction length scale. This can then be solved numerically, and from this solution one can calculate k. 
By changing parameters A and if one can do this for a wide range of velocities, and hen ce obtain a plot of k = k(v). 



2 and a = 3 of our ansatz (3J). It can be seen that a = 3 



This is plotted in Fig. ^ and compared with cases a 
provides an excellent approximation in this regime. 

There is, however, one problem with this simple ansatz, namely that it would give fc(0) = 1. Even though one 
might naively expect this to be the correct limit, it is not so. This can be seen easily as follows. Assume that the 
velocity of say a loop is determined only by curvature, that is, neglect the momentum contribution. (This should be 
valid in the non-relativistic case.) Then k will be approximately given by 



fc ^ <|x 2 |(l-x 2 )> _ 

7j(l — V 2 ) 

This quantity can be easily calculated for the analytic helicoidal solution, yielding 

_ 2^21- 2A 2 /3 



l-A 2 /2 



Hence we find in the small amplitude limit A — » that 



k\\y 



2^2 



0.9. 



(3.6) 



(3.7) 



(3.8) 



We also note in passing that in the relativistic limit A — > 1 this same calculation would give k w 4V2(3/tt) « 0.6 
instead of the true value k = 0, which clearly demonstrates the importance of momentum in the relativistic case. 

The final issue to be considered is the transition between the two regimes. The only reliable way of studying this 
issue is through direct measurement in a string network simulations with ultra-high resolution. We shall report on 
the details of this elsewhere jli| . Here we simply point out that we do confirm the value k nl as the non-relativistic 
limit. It is then easy to find a fitting function for the transition between the regimes which has the correct asymptotic 
limits described previously, that is, 



^ = ^(1-^(1 + 2^)1-1! 



(3.9) 



The additional factors are required to reproduce both the relativistic and non-relativistic limits accurately. Note that 
if one is only interested in the relativistic regime (say for GUT-scale cosmic strings, as in the present paper) then the 
simpler expression (|2.7|), that is 



k Ie i(v) 



2V2 1 - 8v e 



1 + 8v 6 ' 



(3.10) 



should be sufficiently accurate to provide reliable results. On t he o ther hand, a reliable approximation for small 
non-relativistic velocities in the friction-dominated limit extends (3.8) as 



k m (v) 



2V2 



(1~« 2 ) 



(3.11) 



We plot these three ansatze in fig. |5[ and also confirm the validity of fc re i and A: nr in the appropriate limits. 

Before we end this section, however, it is wise to discuss the interpretation of parameter k in this model. It should 
be kept in mind that this is, ab initio, a phenomenological parameter, which accounts for a number of non-trivial 
effects related to the presence of small-scale structures on the strings. By construction, our model does not explicitly 
account for these small-scale effects, and hence they end up somehow encoded in k. One should not therefore infer 
too much from the aesthetic qualities of the function k(v) we find — it is simply a phenomenological parameter that 
does a good job. Presumably this parameter will have a much clearer interpretation in the context of a proper wiggly 
string evolution model p4] , p| . 
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IV. THE EFFECT OF RADIATION BACK-REACTION 



We now turn to some further extensions of the VOS model. The effect of gravitational back-reaction on the long- 
string network jll],[l3| can be included in the evolution equation for the correlation length (2.2) in the same way as 
previously achieved for the evolution of the length of a string loop j!| . For gravitational radiation the following term 
can be added to the right-hand side of ([2. 3D 



dL 
~dt 



8S e 



8TG/j,v 6 



(4.1) 



Here, f is a constant which is a long-string analogue of the T « 65 found for the radiative decay of strings. Of course, 
f will be affected by a number of physical factors such as the presence of small-scale structure, but we can expect 
it to satisfy r < T and it would be surprising if it were very much smaller. Clearly, due to the high velocity power 
(v 6 ) involved in radiative backreaction, this term will not be important in any regime where string motion is strongly 
f rict ion-dominated (and hence non-relativistic) . We note also that there is an interesting coincidence in the ansatz 
(2.7) for the momentum parameter k which too has a v 6 power, but we are unsure as yet whether this has any deeper 
significance. 

For global string radiation into Goldstone bosons or axions, the corresponding radiative decay term at a time t will 

be 



8S a 



8Tv° 



2ir\n(t/S) ' 



(4.2) 



where the logarithmic term arises because of the long-range fields of the global string and S is the string width. For 
cosmological GUT-scale strings, the backreaction term for local strings is TG/i ~ 1CP 4 whereas for global strings it is 
about three orders of magnitude larger. 

Note that the velocity equation has no correction at this order due to the gravitational back-reaction effects. Such 
effects arc already included through the string curvature, which acts as a source for the velocity equation (i.e., the 
1/L term), which will be different in this case. 

Remarkably, the inclusion of the back-reaction term does not affect the existence of a scale-invariant attractor 
solution. However, it does of course influence the quantitative values of the scaling parameters, as well as the timescale 
necessary for this solution to be reached. For example, the inclusion of back-reaction can make the approach to scaling 
much faster. 

If the gravitational back-reaction is non-zero, one can distinguish two asymptotic cases. Firstly, if £ is small (of 
order unity at most) then the effect of back-reaction on the scaling solution will also be small. This will be the case, 
for example, for most local or global string networks in a cosmological context. We can express this as 



7 2 «7o 2 (l + A) 



(1-A) 



(4.3) 



where 70 and vq are the 
form 



unperturbed" scaling values, given by eqns. (2.12), and the back-reaction correction has the 



8j3v$E = 8/3 



m-0) 



1 5/2 



(3{k + 5) 



(4.4) 



On the other hand, for large enough values of X, the back-reaction term will dominate the evolution equation for 
the string length scale L, and the attractor scale-invariant solution has a different form altogether. It is not possible 
to write this solution in closed form, even expressing k implicitly as above. However, it is possible to write it as a 
series. The dominant term and the first correction take the form 



k 

7 ~ 2/3 



8/3E 



k(l-p) 



(1 + A 1 + ...) , 



(4.5) 



HI -(3) 
8/3S 



1/7 



(1-A 1 + ...) 



(4.6) 



with 



7 



Ai 



2 6 / 



7 7 y 



Mi-d)_ 



5/7 



-2/7 



(4.7) 



In fig. ^ we plot the approach to scaling of some relevant string networks in the radiation and matter eras. The 
different timescales for convergence towards the attractor solution are clearly noticeable. Here, we have chosen initial 
conditions that would correspond to somewhat extreme first and second order phase transitions. Also note that 
for the radiation era we have neglected the effect of friction due to particle scattering, in order to reproduce the 
initial conditions often used in numerical simulations of string networks. For each of the cases above, three curves 
are plotted, corresponding to the values E = (no back -reaction) , £ = 1.25 (close to the maximum value that can 
be accurately described by the scaling solution (4.3-4.4)) and £ = 5.5 (beyond which the scaling solution (4.5-4.7) 
becomes accurate). 

For large £, the effects of back-reaction seen in fig. |^ on could be quite dramatic for the string network density, 
but note that they are much less drastic for the string velocities. In particular, we emphasize that gravitational 
back-reaction alone does not slow down a string network to non-relativistic speeds — only a friction-dominated regime 
can achieve this. 

Interestingly, there has been recent work on numerical simulations of global string networks |25f| which explore the 
strong backreaction regime described by (4.5-4.7). These authors report a surprisingly low string density relative to 
the gauged case. For their expanding universe simulations in the realistic case with periodic boundary conditions, 
they find the following radiation and matter era densities respectively, 



Cad = 0.9 ±0.1, Cmat = 0.5 ± 0.1 . 



(4.8) 



These results are perfectly consistent (within the estimated error bars) with our extended VOS model if we adopt a 
back-reaction parameter 

Sax-sim « 3 . (4.9) 

Indeed, this corresponds to the approximate average value for E ax that one would estimate for simulations of this 
resolution. Present limitations on numerical dynamic range give the upper bound \n(t/S) < 5 (at the end of the 
simulation), implying £ a x-sim > 2 throughout. 



It is important to note, however, that the immediate extrapolation of these results (4.8) to a cosmological context 
would be erroneous. For example, for GUT-scale global strings, at the present time we can expect \n(t/S) > 100, 
which implies the app rop r iate backreaction parameter in this case will be E ax < 0.1. Such cosmic global strings are 
firmly in the regime (4.3-4.4) where backreaction effects are small and, in this case, these would reduce the local 
string density by less than 10%. We make a more detailed comparison with |p5l elsewhere ||. 



V. STRING NETWORKS IN GENERAL FRW SPACETIMES 



In this section we discuss the behaviour of our model in more general FRW universes, and in particular in open 

mm- 



universes 



The evolution equations will obviously be affected by the different behaviour of the scale factor, as given by the 
Fricdmann equation 



H 2 = - = H 2 (n R0 a" 4 + r> M0 a' 3 + fl Q o a~ m ) + |A - Ka- 2 . (5.1) 



Note that we are allowing for curvature, and also for an extra fluid whose energy density decays as a~ m . It should 
be kept in mind that the Friedman equation should, in general, contain a contribution for the string density, since it 
is possible that this becomes cosmologically important. 

However, apart from these effects, one must also include an additional correction due to the curvature ]l8p|]. One 
should note that this is essentially the curvature radius of the strings, L, divided by the radius of spatial curvature of 
the universe, 

^ = ^r- (5 - 2) 

Indeed, after a certain amount of algebra, one finds correction terms that are of the form 
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w = l-(l- VL){HLf 



(5.3) 



Note that SI denotes the total density of the universe. For a universe with a critical density, fi = 1, we have w = 1. 
The evolution equation for the correlation length L now takes the form 



2% = 2HL+±&+cv 
at id w z 



(5.4) 



For simplicity wa have also defined a damping length, including both the effects of Hubble damping and friction, 

(5.5) 



with the friction length scale £ / being defined in . 
Similarly, the velocity equation becomes 



dVr, 



dt 



(5.6) 



Note that these are valid for any cosmological scenario^] We do expect the loop chopping efficiency c to be a constant, 
regardless of the cosmological model, since it is supposed to be reflecting a rather deep and fundamental property of 
the evolution of a network. Indeed, we think that whether or not one finds a constant chopping efficiency can in some 
sense be seen as a measure of how accurately the analytic modelling is reproducing the true dynamics of the network. 

We can now re-examine the question of the existence of 'scale invariant' attractor solutions. Again, scaling solutions 
of the form L cx t, L cx H^ 1 or L cx dn , together with Woo = const, will only exist provided one has 

a(t) oc t& , (3 = const. , < (3 < 1 , 



but now we also require 



O = const. 



(5.7) 



(5.8) 



The simplest example of the second condition is of course a flat, JImo = 1 universe, but there are examples of 
cosmological models which have attractors other than = 1 ^(|. In any case, note that there can be additional 
relations between the values of (3 and £1 for specific models. Writing L = jt as before, the scaling solution is now 
given in the implicit form 
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w 



k(k + c) 



w 



,*(!-/?) 



4/3(1 - (3) ' 0(k + c) ' 

where k is (implicitly) the constant value of k(v) for the appropriate value of velocity, and 



2(1-/3) 



(i - n)/3fc(fc + c) 



i 



(1 - n)f3k(k + c) 



1/2 



1 



(5.9) 



(5.10) 



Again, although it may not be immediately obvious, it can be checked numerically that this solution is well-behaved 
for all sensible values of the parameters. If the two conditions above do not hold, then a scaling solution will not exist. 

We should also mention another cosmologically important solution: in an open universe with Q — > 0, a oc i, the 
asymptotic solution is 



L = At (lnt) 



1/2 



Woo = B (In t) 



-1/2 



A = 



D 



2(1 k nr ) 



1/2 



2.13c 1/2 , 



2c 



1/2 



0.21c- 1/2 , 



(5.11) 



(5.12) 



with k nr given by(3.8). Note that this is not a scale-invariant solution, since H = t and dn = tint. In other words, 



by looking at the network one would be able to determine when the curvature-dominated period had started. 



1 There are some additional subtleties involved when discussing the mechanism of loop production in the case of Minkowski 
space string networks, which make it quite different from any cosmological scenario. We shall discuss this important point 
elsewhere ill 611 . 
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VI. DISCUSSION AND CONCLUSIONS 



In this paper we have presented a modified version of the velocity-dependent one- scale (VOS) model @,||J|,§ 
which depends on a single free parameter, the loop chopping efficiency c. We have tested it against the largest and 

1(| , and we find that it provides a good fit to the large-scale scaling 
properties of the string network in both the radiation and matter epochs, as well as in the transition between the two 
eras — we will describe these tests elsewhere Q . These facts and its intrinsic simplicity make this model particularly 
suited for any analytic or semi-analytic study of cosmic strings where one is only interested in the large-scale features 
of the network. 

We have re-analysed some simple evolutionary properties of cosmic string networks in the light of the VOS model 
and corresponding numerical simulations. An important conclusion to note is that any realistic cosmic string network 
is not scaling at any time from just before the epoch of equal matter and radiation through to the present day. This 
is something that must be properly taken into account particularly when discussing string-seeded structure formation 
scenarios with GUT-scale strings. The extended VOS model is also valid when deviations from scaling are even larger 
at late times in a universe which becomes dominated by curvature or a cosmological constant. 

Finally, we considered the effects of radiation back-reaction on the scaling properties of the long string network, and 
we have shown that although the existence (or otherwise) of a scale-invariant attractor solution will not be affected, 
the quantitative scaling properties can be. In some cases, the suppression of string density can be quite dramatic 
(as we saw for small-scale global string simulations), although the string velocities always remain relativistic. For 
the most part, however, the density of a cosmic string network, whether local or global is only affected slightly by 
radiation backreaction effects. 

Despite the many virtues of the VOS model, we are aware, of course, that the small number of available degrees of 
freedom means that this model is unable to provide a proper description of the small-scale properties of the network; 
these are important in a number of cosmological scenarios (and sometimes even crucial). Nevertheless, we believe 
that the phenomenological parameter k does encode some important small-scale structure effects, though clearly a 
more detailed analytic and numerical study is still required. A number of possible approaches to the problem of string 
small-scale structure have been suggested in the literatur e [T^jr^ | , and our own analysis using Carter's elastic string 
model p3] will be discussed in a forthcoming publication [[L4| . 
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FIG. 1. The complete cosmological history of a GUT-scale cosmic string network. The evolution of 7 = L/t is shown in 
the top panel, while the bottom one shows the network's RMS velocity v x . Time is plotted relative to the epoch of equal 
matter and radiation densities; the plots start at the epoch of string formation and end at the present day. At early times, 
the solid curve corresponds to initial conditions typical of a first-order phase transition, while the dashed one corresponds to 
a second-order transition. At late times the solid curve corresponds to a model with J7 m = 1, while the dashed one is for the 
observationally preferred case f2 m = 0.2, = 0.8 (to be discussed below). Note the deviations to the scaling behaviour. 
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FIG. 2. Same as Fig. oj, but for an electroweak-scale cosmic string network. 
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FIG. 3. A close-upof Fig. bl showing the recent evolution of the GUT-scale string network. Also shown is parameter £, a 
measure of the long-string density. Notice the dramatic chages once the cosmological constant dominates. 
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FIG. 5. Comparing our full ansatz for t he m oment um parameter ( p.Sj ) (solid line), and the simpler expressions for the 
relativistic and friction dominated regimes, ( 3.10 ) and ( 3.11 ) (dashed and dotted, respectively), for relativistic (left panel) and 
non-relativistic (right panel) velocities 
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FIG. 6. The effect of gravitational back-reaction on the approach to scaling of a GUT-scale cosmic string network. Left-side 
panels correspond to the radiation era, while right-side ones are for the matter era. All plots have c = 0.23; the back-reaction 
parameter is respectively E = (solid curves), E = 10 (dashed) and E = 50 (dotted). For each case two curves are plotted, 
corresponding to initial conditions typical of a first-order (low density and velocity) or second-order (high density and velocity) 
phase transition. The effect of friction due to particle scattering has been neglected, in order to mimic currently existing 
numerical simulations. 
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